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, Two-dimensional electrons in graphene are known to behave as massless fermions with Dirac- Weyl 

' type linear dispersion near the Dirac crossing points. We have investigated the collective excitations 

' of this system in the presence or absence of an external magnetic field. Unlike in the conventional 

two-dimensional electron system, the v = ^ fractional quantum Hall state in graphene was found 
to be most stable in the n — 1 Landau level. In the zero field case, but in the presence of the spin- 
, orbit interaction, an undamped plasmon mode was found to exist in the gap of the single-particle 

, continuum. 

■ PACS numbers: 
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I. INTRODUCTION 



In recent experimental work, it has been possible to extract an atomically thin, two-dimensional (2D) sheet of 



^ , graphite - graphene, by micromechanical cleavage^. A two-dimensional electron system in graphene exhibits many 

a remarkable properties. In the band structure calculations where electrons are treated as hopping on a hexagonal 
I lattice^, one finds a unique linear (relativistic type) energy dispersion near the corners of the first Brillouin zone 



where the conduction and valence bands meet. As a consequence, the low-energy excitations follow the Dirac- Weyl 
equations for massless relativistic particles^. In an external magnetic field, the spectrum develops into Landau levels, 
S ■ each of which approximately fourfold degenerate^. Recent discovery of the quantum Hall effect in graphene^iii^ has 
' resulted in intense activities^ to unravel the electronic properties of graphene that are distinctly different from the 
, conventional (or, as popularly called the 'non-relativistic') 2D electron systems in semiconductor structures. In this 
Q ' paper, we report on our investigation of the collective excitations of the 2D electron gas in graphene in the presence 
O , or absence of an external magnetic field. 

Graphene has a honeycomb lattice structure of sp^ carbon atoms, with two atoms, A and B per unit cell (i.e., 
I a two-dimensional triangular Bravais lattice with a basis of two atoms). Each atom is tied with its three nearest 
J> . neighbors via strong a bonds that are in the same plane with angles of 120°. The tt orbit {2pz) of each atom is 
' perpendicular to the plane and overlaps with the tt orbitals of the neighboring atoms that results in the delocalized 
TT and TT* bands. There is only one electron in each tt orbit and the Fermi energy is located between the tt and tt* 
^ ■ bands. The separation distance between the nearest neighbor atoms is Ccc = 0.14 nm while the lattice constant is 
^ , a = -v/Sflcc = 0.246 nm. The dynamics of electrons in graphene is described by a nearest-neighbor tight-binding 
\^ • model^ that describes the hopping of electrons between the 2pz carbon orbitals. The first Brillouin zone is hexagonal 
[ and at two of its inequivalent corners (the K and K' points) the conduction and valence bands meet. Graphene is 
■ often described as a two valley (K and K') zero-gap semiconductor. The two valleys correspond to two chiralities of 
' the Weyl-Dirac fermionsiS. Near these two points (the so-called Dirac points), the electrons have a relativistic- like 
^ . dispersion relation, Sk ~ ±?iu|k| and obey the Dirac-Weyl equations for massless fermions. At the vanishing gate 
I ' voltage, the system is half-filled and the Fermi level lies at the Dirac points. 
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] we discuss the properties of the fractional quantum Hall states that reflect the nature of electron correlations in the 
. . 1 system in the presence of a strong magnetic field. We found that the linear dispersion of electrons discussed above, 
leads to a noticeable change in the behavior than what is expected in a conventional two-dimensional electron system. 
In the last section, we discuss the zero-field case using the random-phase approximation and explore the properties 
of plasmons in such a system. 



II. COLLECTIVE MODES 

In the following sections, we describe results of our work on the collective excitations of electrons in graphene. First, 
we discuss the case of the collective modes in the presence of a strong perpendicular magnetic field. In particular. 
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A. In a magnetic field 

An external magnetic field has a significant influence on the energy spectrum of the 2D electron system in graphene. 
Details of the single-electron case has been widely reported in the literature^i^. 



Landau levels 



In the contimmm limit the electron wave function in graphene is a 8-component spinor, ^ s.k.on where s = ±1/2 
is the spin index, k = K, K' is the valley index, and a ~ A, B is the sublattice index. Without the spin-orbit 
interactior>iiii^ii^ the spin degrees of freedom becomes uncoupled from the spatial motion and the Hamiltonian of 
an electron can be described by two 4x4 matrices for each component of the electron spin. If we introduce the 
four-component spinor as {'9s,K,A, "^s.k.b, ^'s,/<'',a, ^^s.k'.b) then in the presence of a magnetic field perpendicular to 
the graphene plane the Hamiltonian matrix has the form 



n = v 



TTx — ITT. 











TTx + JTT. 



V 



(1) 



where n=p + eA/c, p is the two-dimensional momentum, A is the vector potential, and v is the velocity of electrons 
in graphene. It is easy to see from the Hamiltonian matrix ([1]) that the valley index in conserved. The conservation of 
the valley index is easily violated in the graphene systems with a short-range scattering impurity potential or in the 
many-body systems with inter-electron interactions. In both cases the scattering of the electron either by an impurity 
or by another electron introduces the umklapp process and a change of the electron valley index. 

Just as for the non-relativistic system the application of a perpendicular magnetic field to the graphene layer results 
in the Landau quantization. Due to the relativistic nature of electrons in graphene the energy spectrum of the Landau 
levels has a unique form, namely, the energy of the nth Landau level is En ~ sgn{n) y^2ehv^\n\B . This is different 
from the non-relativistic electrons, where the Landau level spectrum is of the harmonic oscillator type, i.e. equidistant. 
En oc n. 

In the ideal graphene system the eigenfunctions of the single-electron Hamiltonian ^ are specified by the Landau 
index, n = 0,±1,±2, . . . and an intra-Landau index m, which depends on the gauge. Each Landau level is fourfold 
degenerate due to the spin and valley degrees of freedom. The corresponding wave functions for an electron in the 
valleys K and K' are described by the vectors 



( sgn(n)zl"l-V|n|-i \ 




V 







(2) 



K' 



Cn 



(3) 



\ sgn(n)il"l / 



where C„ = 1 for n = and C„ = 1/V2 for n ^ 0. The two non-zero terms in K,n ['^ K',n) correspond to occupation 
of the sublattice A (the upper term) and the sublattice B (the lower term). Here 4>n is the Landau wave function 
for a particle with the non-relativistic parabolic dispersion relation in the n-th Landau level. From Eqs. ([2|)-(l3|) it 
is clear that a specific feature of the relativistic dispersion law is the "mixture" of the non-relativistic Landau levels. 
This mixture is present only for n 7^ 0. For n = Q the electron in the valley K or K' occupies only the sublattice A or 
B, respectively. For higher Landau levels the electron in each valley occupies both sublattices, A and B. The wave 
functions in the sublattices A and B are the wave functions of the non-relativistic electrons with different Landau 
level indices, i.e. the relativistic wave function is the mixture of non-relativistic wave functions of different Landau 
levels. As we shall see below, this property of the relativistic electrons strongly modifies the inter-electron interaction 
within a single relativistic Landau level. 
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FIG. 1: Pseudopotentials calculated from Eq. Q are shown as a function of the relative angular momentum for relativistic 
and non-relativistic 2D electrons for the first two Landau levels. The energy is measured in units of ec- 



2. Inter- electron interaction 



In what follows, we shall consider only the partially occupied Landau levels with fractional filling factors. In this 
case the ground state of the system and the excitation spectrum are fully determined by the inter-electron interactions, 
which are completely described by the Haldanc pscudopotcntialsii Vm- Haldanc pseudopotentials are the energies of 
two electrons with relative angular momentum m. The pseudopotentials for the n-th Landau level can be presented 



-.14 



^(/^(q) {FrM)f L„Xq')e 



(4) 



where Lm{x) are the Laguerre polynomials, V{q) = 2TTe'^ / (nlq) is the Coulomb interaction in the momentum space, 
K is the dielectric constant, and Fn{q) is the form factor corresponding to the n-th Landau level. The main difference 
between the relativistic and the non-relativistic electrons is in the expression for the form factor, For relativistic 

electrons the form factor is given by the equation o^^'^^ 



Fo{q) = Lo 



Ln-l 



(5) 
(6) 



while for the non-relativistic particles the form factors in Eq. ^ arc 

F„(q)=i„(gV2). 



Comparing these non-relativistic form factors [Eq. (|II A 2p ] with Eqs. (O-®, we see that the intcr-clectron interactions 
for the relativistic and the non-relativistic electrons arc the same for n — and different for n > 0. 

In Fig. 1 the pseudopotentials calculated from Eq. ^ for the relativistic and the non-relativistic cases are shown. 
For n = the non-relativistic and the relativistic pseudopotentials are the same. The main feature of the pseudopo- 
tentials in the zeroth Landau level, Vm \ is their monotonic decrease with increasing relative angular momentum, 
TO. Comparing Vm^ to the non-relativistic Vm^ we can clearly see that (i) in higher Landau levels the pseudopo- 
tentials become a non- monotonic function of to and (ii) the interaction strength is suppressed at to = and 1, i.e. 

< Vm \ and enhanced at to > 1. i.e. Vm'' > Vm\ From this behavior we can make the predictions about the 
stability and excitation gaps of the FQHE. For example, for the ly = 1/3-FQIIE the main parameter which determines 
the formation the incompressible liquid is the ratio of the pseudopotentials at to = 3 and m = 1, i.e., V3/V1. The 
smaller the ratio the larger the excitation gaps, and consequently the more stable are the FQHE states. Comparing 
the pseudopotentials of the non-relativistic electrons at n = 1 and ?i = we can conclude that the FQHE gaps in the 
first Landau level, n = 1 should be smaller then the corresponding gaps in the zeroth Landau level, n = 0. 

The behavior of the pseudopotentials for the relativistic electrons in the higher Landau levels is similar to that of the 
relativistic ones. The exception is the properties of the pseudopotential in the first Landau level, n = 1. Wc notice in 
Fig. 1 that the pseudopotentials in the first Landau level (n = 1) is larger than the corresponding pseudopotentials in 
the zeroth Landau level, n = 0, for all values of the relative angular momentum, i.e. V^P > Vm'^ . The pseudopotential 
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FIG. 2: Pseudopotentials are shown as a function of the relative angular momentum for relativistic electrons in the first three 
lowest Landau levels. The energy is measured in units of ec- 

in the n = 1 Landau level is also a monotonic function of m. This is different from the non-relativistic case. Another 
important property of the pscudopotcntial of the relativistic system is that the ratio V3/V1 is the smallest in the first 
Landau level, n ~ 1. In Fig. 2 we present the results for the pseudopotentials of the relativistic system in the lowest 
Landau levels. From this figure we can also see that the ratio V3/V1 is the smallest for n = 1. Based on this property 
of V^ii"'' we can conclude that the largest gap of the 1/3-FQHE state should be expected in the first Landau level, 
n = 1. Therefore the FQHE in the relativistic system should be easier to observe in the first Landau level, but not in 
the zeroth Landau level as in the case of the non-relativistic electrons. 

Due to the presence of both the spin and the valley degeneracies the graphene system becomes more complicated 
than the standard non-relativistic electron system, where only the spin degree of freedom is present. In a magnetic 
field the spin degeneracy can be lifted due to the Zeeman energy, leaving the Landau levels of the graphene system 
doubly degenerate due to the valley index. Mathematically, the system then becomes equivalent to a double layer 
non-relativistic system, where the valley index is the layer index. Usually, in the double-layer system there is an 
asymmetry in the interaction Hamiltonian. This means that the interaction strength between the electrons in the 
same layers is different from that in different layers. This is due to a finite separation between the layers. Comparing 
this property of a double layer system to the graphene system we can say that the graphene system is a double-layer 
system with zero separation between the layers. So the system is completely SU(2)-pseudospin symmetric, where the 
pseudospin is associated with the valley index. But this is not entirely true. There is an asymmetry in the graphene 
system as well. One of the types of asymmetry is related to the lattice structure of graphene — the interaction between 
the electrons in the same sublattice is stronger than the interaction between the electrons in different sublatticcs^. 

To estimate the corrections due to the asymmetry related to the lattice structure of graphene we assume that 
the electrons are localized at the discrete points corresponding to the lattice structure of the graphene. Then the 
interaction matrix element, ^i,j"i,i2,j2J between the single-electrons states, ipi, can be expressed in the following form 

V^,,n,^2,J. = (ri)V*,(r2)l^(ri - r^)^/.,, (rs)^'.^ (ri) (7) 

ri r2 

where the index i in the wavefunction ^pi indicates collectively the Landau level, valley, sublattice, and intra-Landau 
level indices. The sums in (O run over the discrete positions of the electrons corresponding to one of the sublatticcs 
of graphene. In the continuum limit the sums should be replaced by the integral. Since two sublattices of graphene 
are shifted by the vector Tq = a(0, l/>/3), in the continuum limit the expression for the interaction matrix element 
between the states belonging to different sublattices should contain not the ^(ri — but ^(ri — r2 — Tq). This 
introduces the difference in the interaction strength between the electrons in the same sublattice and in different 
sublattices. In terms of the pseudopotentials this means that if the pseudopotential is calculated between the states of 
different sublattices then V{q) in the expression ^ should contain an additional factor exp(iq • Tq) and the integral 
over q should be replaced by the 2D integral over q. 

In the zeroth Landau level, n = 0, the electrons in the valley K occupy sublattice A only, while the electrons in 
the valley K' occupy sublattice B. Then the pseudopotentials corresponding to the interaction between the electrons 
belonging to the same valley are given by the expression ^ without any modifications. The expression for the 
pseudopotentials corresponding to the interaction between the electrons in the different valleys should contain an 
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additional factor and can be written as 



(0) 



K,K'. 



V^ 



(0) 

K.K,-n 



y(g)e*i-"-«/'i™((z2)e-9^ 
q^dq 



47r 



V{q)L„,iq')e-'^ . 



(8) 



Since tq = the asymmetric correction is proportional to a small parameter (a/l)^. 

In the higher Landau levels the electrons in the valleys K and K' occupy both sublattices A and B. This results in 
the form factor of ■|(L„ + in Eq. (jl]). The coefScients L„ and L„_i belong to different sublattices. This can 

be schematically presented in the following form: (i) for the intravalley interaction, we have 



i [L„(A)i„(A) +L„_i(B)i„_i(B) + 2L„(A)L„_i(B)] 



and (ii) for the intcrvalley interaction 

1 



[L„(A)i„(S) + L„_i(B)L„_i(A) + 2L„(A)i„_i(A)] 



(9) 



(10) 



Then following the same procedure as for the n = Landau level we obtain the expressions for the pseudopotentials 
in the higher Landau levels as 



in) 



K,K,m 



yin) 



^,V[q)L^[q-)e-^' 



and 



V^ 



(") 



K,K'm 



J ^Viq)LUq')e-'^' F^iq)e^^-'^°^^ + ii„L„_i (l - e'^''^^^ 



It is convenient to rewrite the expressions pT|) and as 



^K.K' 



(n) 



K.K,i 



(4^ 



V{q)L„,{q^)e~''" [L„ - L„_i]' (l - e'l-"-"/') . 



(11) 



(12) 



(13) 



Similar to the zeroth Landau level the asymmetry correction in Eq. (I13|) is proportional to a small parameter (a/l)'^. 

Since the Coulomb interaction between the electrons in graphene does not conserve the valley index there is an other 
mechanism for violation of the SU(2) valley symmetry of the graphene system. This mechanism is related to the inter- 
valley scattering. In the lowest order in the small parameter a/ 1 the main scattering process is the backscatteringi^. 
During this process the electron from the K valley is scattered into the K' valley, while the electron from the K' 
valley is scattered into the K valley. The interaction matrix elements corresponding to this process arc determined 
by the pseudopotentials ^ where the Coulomb interaction in the integral should be replaced by ^(q + /AK). Here 
AK = K — K' = (27r/a)(— 1/3, l/VS). Therefore the pseudopotentials describing the backscattering process are given 
by the expression 



Vj,t^ J ^Viq + IAK) [F„(g)]^L„(g2)e-^^ 



(14) 



The leading order term in V^"jj can be found by simply replacing y(q + / AK) by F(/AK) oc {a/ 1). 



3. FQHE in graphene 

With the pseudopotentials for Dirac electrons at hand, we now evaluate numerically the energy spectra of the 
many-electron states at the fractional fillings of the Landau level. The calculations have been done in the spherical 
geometryi^ with the pseudopotentials given by Eq. In the spherical geometry the radius of the sphere R is related 
to 25* of magnetic fluxes through the sphere in units of the flux quanta as i? = ^fSl. Here 2S is an integer number. The 
single-electron states are characterized by the angular momentum S, and its z component, Sz- Therefore, at a given 
magnetic field, i.e for a given flux 25*, the number of available states in a sphere is (25* + 1). Then for a given number 
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FIG. 3: (a) The energy spectra of an eight-electron v = 1/3-FQHE system shown for different Landau levels: n = (stars) 
and n = 1 (filled circles). The flux quanta is 25 — 21. (b) Energy spectra of the six-electron v = 1/5-FQIfE system is shown 
for different Landau levels: n — (stars) and n = 1 (filled circles). The flux quanta here is 2S = 25. 



of electrons N the parameter S determines the fiUing factor of the Landau level. Due to the spherical symmetry of 
the problem, the many-particle states are described by the total angular momentum L and its z component, while 
the energy depends only on L. The energy spectra of a many-particle system is found by the standard procedure of 
calculating numerically the lowest eigenvalues and eigenvectors of the interaction Hamiltonian matrisi^. 

It was shown in the previous section that based on the analysis of the Haldane pscudopotcntials we can conclude 
that the FQHE states in graphene should have the largest gaps in the n = 1 Landau level. Here we check this 
statement for ly = 1/m incompressible states. In the spherical geometry, such states are realized at S' = {m/2){N — 1). 
The 1 /m state in the higher n-th Landau level is defined as a state corresponding to the 1 /m filling factor of the n-th 
Landau level, while all the lower energy Landau levels are completely occupied. If the electron system is fully spin 
and valley polarized then we should expect that the ground state to be the Laughlin statei2i2£ which is separated 
from the excited states by a finite gap. 

In Fig. 3(a) the calculated energy spectra are shown for the 1/3-FQHE state and for different Landau levels. Since 

the relativistic pseudopotentials, Vm \ for n = Landau level is similar to that of non-relativistic case, the 1/3 state 
and the corresponding energy gap will be the same in both cases. The deviation from the non-relativistic system 
occurs at the higher Landau levels. From Fig. 3(a) we can clearly see that the energy gap of the 1/3-state at the n = 1 
Landau level is enhanced when compared to that of the n = Landau level. At higher Landau levels, i.e., for n > 1, 
the excitation gaps are suppressed, which means that at the n = 1 Landau level the electron system in graphene 
has the strongest interaction with the largest incompressible gapiSi. This is different from the non-relativistic case, 
where the energy gap monotonically decreases with increasing Landau level indexes. At other filling factors of the 
type 1/ ~ 1/m we should also expect the same increase of the energy gaps at the n = 1 Landau levels. The effect 
is however not as pronounced as for v = ^ since at smaller filling factors the pscudopotcntials with a larger relative 
angular momentum, m, becomes important, for which the difference between the pseudopotentials at n = and n = \ 
Landau levels is small [see Fig. 2]. This behavior is illustrated in Fig. 3(b), where the results for the v ~ 1/5 state 
are shown. We can see that the difference between the excitation spectra at n = and n = 1 Landau levels is much 
smaller for this filling factor. 

The results shown in Fig. 3 correspond to the fully spin and valley polarized systems. This polarization is achieved 
at a high magnetic field due to the Zeeman splitting and the valley asymmetry. It is well known that even without 
the Zeeman energy the ground state of the non-relativistic system at = 1/m is fully spin-polarize d^°i^^ with the 
spin equal to 5 = N /2^ where N is the number of particles. The situation is the same in graphene: the ground state 
of the V = 1/m liquid is fully spin and valley polarized. The results shown Fig. 3 describe the polarized excitations 
in such a system. Another type of neutral excitations of the incompressible liquid is the spin or pseudo-spin (valley) 
excitations. To explore these excitations we present the results of our calculations of the energy spectra for the z/ = 1/3 
'unpolarized' system. Here we assume that the system is fully spin-polarized but valley-unpolarized. To characterize 
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FIG. 4: The energy spectra of a six-electron valley-unpolarized v — 1/3-FQHE system shown for different Landau levels: n = 
(stars) and n = 1 (dots). The flux quanta is 25* = f5. The spin- wave excitations are illustrated by solid (n = 1) and dashed 
(n = 0) lines. 



the states in such a system we consider the valley index as a pseudospin, r. Results of these calculations are shown 
in Fig. 4 for the zeroth and for the first Landau level. The ground state is at zero angular momentum and has a 
pseudospin equal to r = N/2, i.e. the pseudospin-polarized ground state. One type of excitations in a such system is 
the spin-waves. They are marked by the solid (n = 1 Landau level) and dashed (n = Landau level) lines in Fig. 4. 
Another type of excitations corresponds to the energy branch formed by the lowest states at each angular momentum. 
The specific feature of these states is that the angular momentum of the state is related to the pseudospin value by 
the expression {N/2) — t = L. The physical meaning of these states is the Bose condensation of L noninteracting 
pseudospin waves^^. Clearly, in all the cases the energy scale at n = 1 is larger than at n = 0, which again illustrates 
the stronger interaction effects at n = 1. However, this is not the general rule for the n = 1 Landau level since we 
have seen from Fig. 2 that the pseudopotential at the zero relative momentum is stronger at the zeroth Landau level. 
This pseudopotential becomes important only for the unpolarized states when two electrons occupy the same spatial 
point, i.e. they have zero relative angular momentum. For those states the interaction effects can be stronger at the 
zeroth Landau level, as is the case of the = | incompressible state^i. The excitation gap of the unpolarized v = ^ 
state at the n = Landau level is larger than the excitation gap at the n = 1 Landau level^i. 

In the above calculations we have also included the valley asymmetry terms discussed in the previous section [see 
Eqs. ([8|)- p4|) ]. These corrections result into a very small shift of the energy levels up to magnetic field of 50 Tesla. 
Indeed, these terms are of the order of a/l. The magnetic length at = 50 T is about 3.6 nm, while the lattice 
constant of graphene is a = 0.246 nm. The parameter a/l is then really small (~0.07). 

Finally, from the results presented above we conclude that the 1/m-FQHE state in graphene is most stable at 
the n = \ Landau level. The inter-electron interaction effects are therefore more pronounced at the n = \ Landau 
level^. This tendency is just the opposite to that of the non-relativistic system, where the excitation gap decreases 
monotonically with increasing Landau level index. The enhancement of the interaction effects at the higher Landau 
level index however depends on the filling factor. 



B. Zero magnetic field 



In the pseudospin space, the zero-magnetic-field Hamiltonian of a spin-up electron with a wavevector around the 
K point isiiii^iH Ji ~ vp ■ cr + AgoCz with cr = {ax, <Jy,crz) the Pauli matrices and p the momentum operator. Here 
Ago is the strength of the spin-orbit interaction (SOI). The eigenstates of the Schrodingcr equation iJ^' = E'^i are 

readily obtained as '^/^{r) = e''=-^ -e'^^^'cl^{ak + \i%) ) "^'^^^ '^'^^^'^^ " X./Al^TW^ for A = +1 denoting 

the conduction band and A = —1 the valance band. Here tan0fc = ky/kx,ta,nak = hvk/ Ago, and k = \J kx + ^y- 

Using the techniques developed for the multicomponent systems^ii^, it is straightforward to show that the RPA 
Coulomb interaction in the Fourier space U{q,uj) obeys the equation 

U{q,u}) =vo + votlo{q,u})U{q,u}) (15) 
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FIG. 5: Diagrammatic illustration of the RPA dressed Coulomb interaction. 
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FIG. 6: Plasmon spectrum (thick curves) of an electron gas in an intrinsic graphene {Ef ~ 0) at temperatures T = 1.42 K 
with Ago = 0.08 meV. Intra- (dark shaded) and inter- (light shaded) band single-particle continuums are also shown, ujl and 
ujh are the lower and upper borders separating the white (EHC gap) and shaded areas respectively. 



with the electron-hole propagator 



no('7,w) 



E 

A.A'.fc 



I A. A' 



(9) 



fiE'k+J-flE^] 



Ek+q 



E'k 



i5' 



(16) 



as illustrated by the Feymann diagram in Fig. [S] Here wq = /{2eQ€iq) is the two-dimensional Coulomb interaction 
(in Fourier space) with the high-frequency dielectric constanl^ = 1 and g^'^ (q) is the interaction vertex. 

The factor four in Eq. (|16p comes from the degenerate two spins and two valleys at K and K'\ the vertex factor reads 

\9k^ (9)1^ = [1 + cosafc+q cosafc -I- AA' sinafc+q shiak{k + qcos9)/\k + q\]/2 with 9 being the angle between k and 
q. Since the chiral property of the system prohibits the intra-band backward scattering at q = 2k and the inter-band 
vertical transition at q = under the Coulomb interaction in the system, we have |gfc'~'*'(0)p = |g^''^(2fc)p = 0. 
The collective excitation spectrum is obtained by finding the zeros of the real part of the dielectric function e{q, uj) = 
1 - Vo{q)tlo{q,uj). 

In the presence of the SOI, an energy gap opens between the conduction and valence bands and the semimetal 
electronic system in graphene is converted into a narrow gap semiconductor system. At the same time, a gap is 
opened between its intraband single-particle continuum uj < lul = hvq and its interband single-particle continuum 
^ > = 2^ -|- h?v'^q'^ /A. However, the system differs from a normal narrow gap semiconductor due to its 
peculiar chiral property. In this paper, we have chosen the magnitude of the SOI strength to be around 0.08 — 0.1 
meV in grapheneiii^. The result can be easily applied to Dirac gases with different Ago by scaling the energy and 
wavevector in units of Ago and fcgo = Aso/(?iw) respectively. 

At zero temperature or for T <C Ago, the intraband transition is negligible and > 0. There is no plasmon mode 
in the system. With an increase of the temperature, holes appear in the valence band and electrons in the conduction 
band. The intraband transitions are enhanced and contribute to the electron-hole propagator of Eq. (jl6[) and a dip in 
Cr at the intra-band EHC edge lol. This dip in results in plasmon modes above uj^. For Ago = where loh = lj^, 
the intraband (interband) single-particle continuum occupies the lower (upper) part of the lj ^ q space below oj^ 
(above ljl) and the plasmon mode are Landau damped. In the presence of the SOI, i.e. for Ago 7^ 0, a gap of width 
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FIG. 7: Temperature dependence of the real part of the dielectric function at the edges of the intra- and intersubband single- 
particle continuum ujl (dotted curve) and u>h (solid curve) at g = 0.05 x 10^ cm~^. 



ljh — ojl is opened between the intra- and interband single-particle continuum and an undamped plasmon can exist 
in this gap as shown in Fig. [S] This plasmon mode may perhaps be observed in experiments. 

The appearance of the undamped plasmon mode in the presence of the SOI is a result of the interplay between the 
intra- and the inter-band correlations which can be adjusted by varying the temperature of the system in experiments. 
To show the temperature range in which an undamped plasmon mode exists, in Fig. [7] we plot Cri^L) (dotted curve) 
and er(wij) (solid curve) as functions of the temperature T aX q = 0.05 x 10^ cm~^. For Ago = 0.08 meV, an increase 
of the temperature from T = leads to an increase of the ratio of the intra- to the inter-band correlation while in 
the EHC gap (w^, ^ w ^ uih) decreases and crosses zero. There is no undamped plasmon mode when the inter-band 
correlation dominates at T < 1.1 K and when the intra-band correlation dominates at T > 3.3 K. In the temperature 
regime 1.1 K < T < 3.3 K or T w 2Aso when the intra- and inter-band correlations match, however, eri^h) < while 
er.(wjy) > and one undamped plasmon mode exists. 

In summary: calculating the dynamic dielectric function taking into account the intrinsic spin-orbit interaction in 
graphene, we have studied the collective excitations in graphcne. The Dirac electronic system in graphene is converted 
into a narrow gap semiconductor with chiral property by the spin-orbit interaction. As a result, an undamped collective 
excitation was found to exist in the spectral gap of the single-particle continuum and is perhaps observable in the 
experiments. More detailed results can be found elsewherei^. There have been a steady flow of reports in the 
literature on the electronic properties of graphene. Interestingly, our SOI-dependent dielectric function has recently 
been employed to explore the possibility of Wigner crystallization in graphene^. 
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